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ABSTRACT 

We derive general equations for non-linearity corrections and statistical uncertainty (variance) 
estimates for data acquired with near-infrared detectors employing correlated double sampling, 
multiple correlated double sampling (Fowler sampling) and uniformly-spaced continuous readout 
techniques. We compare our equation for the variance on each pixel associated with Fowler 
sampling with measurements obtained from data taken with the array installed in the near- 
infrared cross-dispersed spectrograph (SpeX) at the NASA Infrared Telescope Facility and find 
that it provides an accurate representation of the empirical results. This comparison also reveals 
that the read noise associated with a single readout of the SpeX array increases with the number 
of non-destructive reads, n r , as n" 16 . This implies that the effective read noise of a stored 
image decreases as n~ a34 , shallower than the expected rate of n~°' 5 . The cause of this read 
noise behavior is uncertain, but may be due to heating of the array as a result of the multiple 
read outs. Such behavior may be generic to arrays that employ correlated or multiple correlated 
double sampling readouts. 

Subject headings: instrumentation: detectors — methods: data analysis 

1. Introduction 

Correlated double sampling (CDS) and multiple correlated double sampling (MCDS, also known as 
Fowler sampling; Fowler & Gatley 1990) are common techniques used to read out near-infrared arrays. 
These techniques are used, for example, for recording data obtained with the arrays in SpeX at the NASA 
Infrared Telescope Facility (Rayner et al. 2003), in NIRSPEC at the Keck telescope (McLean et al. 1998), 
in the IRCS at the Subaru Telescope (Kobayashi et al. 2000), and in ISAAC at the Very Large Telescope 
(Moorwood 1997). An alternative readout technique, known as continuous, or "up-the-ramp" sampling, was 
used for recording data obtained with the CGS2 spectrograph on the United Kingdom Infrared Telescope 
(Chapman et al. 1990). A variant of this technique is the recommended readout scheme for the detectors in 
the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) aboard the Hubble Space Telescope 
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(Roye et al. 2003). The relative advantages and disadvantages of these techniques in various noise- limited 
operating regimes have been summarized by Garnett & Forrest (1993). 

As part of the effort to develop a set of semi-automated routines to reduce SpcX data (Spextool; see 
Cushing et al. 2003), we have derived equations for the non-linearity corrections and the uncertainties on 
each pixel value recorded by an array that employs CDS, MCDS, and uniformly-spaced continuous sampling 
readout techniques. In this paper, we present a simple model for the CDS and MCDS readout procedures of 
a near-infrared array, define the relevant parameters and derive equations for the non-linearity corrections 
and variances on pixel values associated with these techniques (§2). The corresponding equations for the 
uniformly-spaced continuous sampling technique are presented in §3. Verification of the MCDS equations 
and their application to the SpeX array for the determination of the gain and read noise are presented in 
§4. A summary of our results is given in §5. 

2. Correlated and Multiple Correlated Double Sampling 
2.1. Ideal Detector 

Figure la shows a schematic of the signal recorded on an integrating node as a function of time for an 
ideal detector exposed to a constant source flux and read out using the MCDS technique. The ideal detector 
is perfectly linear over its entire well-depth, thermally stable, and not subject to any time dependent effects 
which vary during integrations or readouts (such as bias drift, image persistence that decays over time, 
or amplifier glow associated with readouts). Excellent descriptions of such systematic effects, which will 
assumed to be either negligible or correctible (and therefore ignored) in our analysis, are given by Skinner 
& Bergeron (1997) and Royc et al. (2003) for the NICMOS array. 

At time t = 0, the array is reset to the bias level b, which is assumed to be constant in time during the 
integrations and readouts. 1 In near-infrared arrays, the bias is a negative voltage set by the electronics and 
therefore is "noise-less" . Because it is assumed to be simply an additive constant that defines a zero point 
for count levels at each pixel, it can be removed from all equations without any loss of generality (as long as 
it does not vary) . We assume it takes a time t reset to reset an individual pixel. The reset procedure may be 
global (all pixels in the array reset at once, as in SpeX) or sequential (array reset occurring pixel-by-pixel, as 
in NIRSPEC). After reset, the array is immediately read out non-destructively n r times. The time needed 
to read an individual pixel is t read , and a readout occurs sequentially. Reading the entire array, therefore, 
requires a total time 5t, given by the time needed to read out a single pixel times the number of pixels read 
rip, plus any overhead or waiting time between reads t wa u, 

fit ^ptread "i" ^wait • (1) 

In SpeX, for example, the 1024 x 1024 array is divided into quadrants, which are read simultaneously and 
each of which has 8 readout outputs. For this array, t reset ~ 30^s, t rea d 10^s, and n p — 2 15 ; therefore fit 
is on the order of 300 ms. 

The first set of reads provides an estimate of the mean pedestal level p. After an exposure time At, 
the array is again read out n r times, which provides an estimate of the mean signal level s. (In the CDS 



x Variables written in boldface denote two-dimensional quantities, stored in images, which may have different values at each 
of the pixels in the image; operations performed on images are carried out on all pixels individually. The individual values at 
a given pixel k, I arc written in italics with subscripts, e.g., s^;. 
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technique, n r = 1.) Therefore, At is given by the time difference between the start of the first read of the 
pedestal pi and the start of the first read of the signal si. In most implementations of these techniques, the 
values of neither the individual readouts nor the mean pedestal and signal levels are stored; instead, usually 
only the net source counts (in DN) for the n r non-destructive readouts over an integration time of At, given 

by, 

Snet = Sl ~ Pi ( 2 ) 

i=l »=1 

= n r s — n r p (3) 

is recorded. The estimated source flux, I (which includes contributions from the dark current, background, 
sky, and object), in units of DN s _1 is then given by, 

'-Mi*-^ ' (4) 

\i— 1 i—1 / 

Hence, I is the slope of the solid diagonal line in Figure la. 

The goal is to determine the uncertainty, or variance, in the estimated source flux I at each pixel. As 
can be seen from Equation 4, and will be further demonstrated below (§2.4), expressions for the individual 
pedestal and signal readout values are needed for the computation of the variances in I. From Figure la, it 
can be seen that for a linear detector, after a reset the counts (in DN) recorded at the end of the ith read 
of the pedestal, Pi, and the signal, s,, can be expressed in terms of the recorded net source counts S„ et or 
the estimated source flux I as 



,. ~Snet St 

0- f (5) 

n r Lxt 

(i-f)ISt . (6) 



and 



s, = ^*+(i-f)^*^ (7) 
n r n r At 

= + (8) 

n r 

= IAt + Pl , (9) 

respectively. Here f is a fraction, f < 1, whose value at a given pixel depends on the position of the pixel 
in the readout sequence. The quantity (i — fk,i)St is the time after the reset when the pixel k, I is read out 
for the zth time. If a total number of n p pixels are read out, / for the mth pixel in the readout sequence is 
given by 

j» -. t rea d t rese t 

f = 1 - m —-— (10) 

for a global reset and 

p tread , . A preset /1 1 \ 

for a sequential reset. 
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A common observing technique employed with near-infrared arrays is to immediately repeat an integra- 
tion n c times and "co-add" the results. In this case, the entire array addressing procedure (reset, n r pedestal 
reads, n r signal reads) is repeated n c times. Usually, the net source counts resulting from the individual 
integrations are not recorded. Rather, the total source counts given by 

Stat = ^ S ne t,j = Si j - ^ ^ Pi j (12) 

j=l j=l i=l j=l i=l 

= n c n r s — n c n r p . (13) 
is stored. The estimated source flux I in units of DN s _1 is given by, 

I = — ^rrS tcrf . (14) 
n c n r At 



Therefore, 



and 



,. p-, Stot St n r\ 

P ' = (j - f) ^Al ' (15) 



Stot . ,. #»\ Stot St 

Si = + (*-f) -TT (16) 

n c n r n c n r At 



(17) 



2.2. Non-Linear Detector 

Unfortunately few astronomical detectors are perfectly linear over their entire well depth. Near-infrared 
arrays are inherently non-linear devices because the detector capacitance varies during an integration as 
photons are detected, (see e.g., McCaughrean 1988). Figure lb shows a schematic of the signal recorded 
on an integrating node as a function of time for a non-linear near-infrared array read out using the MCDS 
technique. The recorded signal (the solid line) falls below that expected in the ideal case (the diagonal 
dashed line). In the absence of other systematic and time variable effects (such as image persistence, bias 
drift, or thermal instability of the array), the ratio of the theoretically-expected signal to the recorded signal 
yields the non-linearity curve of the array, which increases steadily from unity as a function of recorded 
counts, 

C nl ( Sl ) = ^ , (18) 

where s\ m is the signal recorded by an ideal (linear) detector. Therefore, the values of the individual readouts 
of the pedestal and signal are given by 

s 2 = -^— and Pi = ^L . (19) 

If the non-linearity of the detector is not accounted for, the estimated source flux I can be substantially less 
than the true source flux. Furthermore, as can be seen from Figure lb, the values of the individual pedestals 
and signal readouts derived from Equations 5-9 will be underestimated. 
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The corrections C„; for the non-linearity of optical and near-infrared arrays are often derived from flat- 
field data acquired expressly for this purpose. A series of flat field images taken with gradually increasing 
exposure times are used to generate a curve of the recorded signal as a function of the exposure time, or 
equivalently estimated count rate as a function of recorded counts. (Of course, all systematic effects, such as 
those mentioned above, must be corrected for and removed before the non-linearity curve can be determined 
from the flat-field data.) The curve can then be fit with a continuous function (e.g., a polynomial), and 
the deviation of the function (at large count values for long exposures) from a straight line (fit to the count 
values recorded during short exposures) provides a measure of the non-linearity correction C n i(sk,i) as a 
function of the detected counts Sk,i in pixel k, I. For the SpeX array, for example, it was found that the 
non-linearity curve can be approximated extremely well at each pixel over the range to ~ 8500 DN (which 
represents ~ 90% of the well depth and is close to saturation) by a rational function, 



C„,(s)« — — 3 , (20) 

1 + ai • x + a 2 • + a 3 • x d 

where the ai, a 2 , and a 3 are constants, x = s — pfiat, and p/; a t is the pedestal level estimated for the flat 
field images (using Equation 23 below). This curve can then be used to correct the signals recorded by the 
detector (after removal of the bias, dark current, and any systematic effects) to the "true" values that would 
be recorded by a perfectly linear detector. 



2.3. Non-Linearity Corrections 

When attempting to apply a non-linearity correction to data acquired with CDS or MCDS techniques, 
it is important to keep in mind that the corrections must be made to both the signal and pedestal values 
separately, not to the net source counts, S ne t or Stot- However, because these readout procedures do 
not record separate signal and pedestal values (see §2.1), we must estimate these a posteriori from S tot . 
Fortunately, this can be done relatively easily, with the definitions given in §2.1, although in general an 
iterative procedure is required. 

We assume that the non-linearity corrections apply to the mean signal and pedestal values above the 
bias level, that is, the relative values above the bias rather than the absolute values. Furthermore, we assume 
the bias level remains constant during the integration and readouts, the dark current has been subtracted, 
and any other systematic effects have been corrected for before the non-linearity corrections are applied. In 
this case, first estimates of the pedestal and signal values can be calculated from 

j=i i=i 

1 ED-^S (22) 



n c n r ^ ^ n r At 

j = l V- 1 



'tot 



St 



n c n r At 



n r + 1 _ 



(23) 



and 



respectively. 



s (1) = J-S^ + pW , (24) 



- 6 - 



These values are then corrected for non-linearity using the C„; curve 

s( 2 )=s( 1 ).C„ ; (s (1) ), P (2) =P (1) .C n! (pW) (25) 

and a new estimate of the net source counts is computed from 

S<£{ = n c n r s (2) - n c n r p (2) . (26) 

It should be noted that the corrected values of s, p, and S tot depend on the number of non-destructive reads, 
n r (see Eq. 23). 

(2) 

Using the corrected estimate of the total source counts S). ot , a new estimate of the mean pedestal level, 
p( 3 ) can be computed from Equation 23. Adding this to the original values of S tot using Equation 24 yields 
a new estimate of the signal, . The new estimates of s and p can then be corrected for non-linearity 

(3) 

using the C n i curve and then used to generate another estimate of true source counts S tot . The procedure 
can be repeated as many times as necessary to achieve a desired convergence criterion. 

For SpeX data we have found that convergence is rapidly achieved and that, unless the mean detected 
signal s is close to the top of the well (> 8500 DN, highly non-linear), the correction to the estimate of 
the true source counts is less than 1% after three iterations; therefore, s=s ^tot- An example of both 
the necessity of making corrections for non-linearity and the effectiveness of the procedure outlined above 
can be seen in Figure 2, taken from the paper by Cushing et al. (2003). This figure shows that spectra 
in adjacent spectral orders, reduced without incorporating the non-linearity corrections, do not match one 
another in either intensity level or spectral slope in the overlapping wavelength region. After the non-linearity 
corrections described above are made, however, the spectra agree in both intensity and slope to better than 
a few percent over the entire overlap region. 



2.4. Variance Estimates 

Because the array is read out non-destructively in the CDS and MCDS (as well as the "up-the-ramp" ) 
techniques, the values recorded at a given pixel for each readout are correlated. Therefore, the variance V/ 
of the estimated source flux I cannot be estimated simply by the photon noise and the read noise for the 
individual readouts added in quadrature. The covariance between any two reads must be included in the 
variance estimate (Garnett & Forrest 1993). 

We start by considering the case where n c =l. The variance of I, V/, is given by the standard error 
propagation formula (Bcvington 1969; Garnett & Forrest 1993), 



E 

i=l L 



31 

dsi 



dl 

1 \dpi 



(27) 



j=2 i<j 

+2EE 

i=i j=i 



cov( Si , Sj )^ (|1 



cov(pi, Pj ) [ — 



dl 
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where cr 2 . is the variance of the zth readout value and cov(xi,xj) is the covariance between the ith and jth 
read. 

The variance of the ith readout value is given by the sum of the square of the photon noise and the 
square of the read noise. The read noise, er read} is usually specified in terms of electrons. If the array has 
a gain g (with units of electrons per data number, e~ DN -1 ), conversion to DN gives a variance associated 
with the read noise of f 2 ead /g 2 - F° r Poisson noise, the variance of each readout value, expressed in DN 2 , is 
given by 

a 2 s = * and al = E* . (28) 
g g 



The covariance between the zth and jth readout is given by the square of the photon noise of the zth 
read for j > i (Garnett & Forrest 1993). This can be easily shown as follows. For any two reads z, j, with 
j > i, the associated readout values are x, and Xj, which are related by 

xj = x, + Aj_i (29) 

where Aj-i is the difference in counts between the two reads. Then 

COv{*i,Xj) = ((x; - (x;)) • ( Xj - ( Xj ))) (30) 

= (x 2 ) - (x 4 ) 2 + (xi • Aj_i) - (xi) • (Aj-i) 
= cov(xi,Xi) + cov(xi, Aj_i) 

= < (31) 
where the angle brackets ( ) denote the average. Therefore, we have 

cov(si,Sj) = (T Sz and cov(pi,pj) = a p Jor j>i (32) 

and 

cov(Pi,Sj) = a Pi . (33) 
Upon substitution, Equation 27 becomes, 



V, = 



' j + +2 wK Hi _ 2 yyPi (34) 

t1\s g 2 g g 2 / f^>t1\s g/ ~1t^ s I 



n?At 2 



j=2 i=l x ° ° 7 i=l j = l 



1 " r ™ r ™ r 2n I 

, ^(s 4 +p l ) + 2^(n r -i)(s l +p l )-2n r .^ Pl + — ^cr 2 eQ(; l (35) 



1 x x -> , \ ^ 2n r 



^(s i + p0 + 2n r ^s i -2^z(s i + p i ) + — -<r 2 ead [ . (36) 



grz 2 Ai 2 ^ ^ g 



We now assume that the array has a perfectly linear response, such that Si = s- m , pi = p\ m , and 
S ne t = S^" t . If the array is not linear, the signal and pedestal images must be first corrected for non- 
linearity as outlined above. Explicitly accounting for non-linearity would introduce factors of C„i(sj) and 
C„;(pi) in the equations given below, which would result in an increase in the variance of any given pixel value 
beyond that expected for an intrinsically linear array and would also prevent us from deriving completely 
general equations for the variance in terms of recorded quantities. However, as we will demonstrate below 
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(§4) , this increase can be easily accounted for by modifying the values of the gain and read noise from their 
"intrinsic" values. 

Substituting the definitions of Sj and pi for a linear, stable (constant b) detector from Equations 5 and 
7 into Equation 36, we have, 



1 



gn 2 At 2 

- 2 E* 



E 



n r 



+ (i-f) 



2Sr,, R + St 



n r At 



2n r 



Snei ,^ £\ ^net St 

n r n r At 



i=l 



+ (i-f) 



2S ne j <5t 
n r At 



2n 



— /T 2 

^ read 



Using the relations 



n _^ n _^ 

^i=-n(n+l), ^ 2 = -n(n + l)(2n+l) , 



8 = 1 



we can simplify Equation 37 as follows, 



'net 



gn r At 2 
I 



1 <5t (n 2 - 1) 



gAt 



3 At n, 

2 



+ 



o_2 

ZCr read 

g 2 n r At 2 



. 1 ft (n 2 - 1 
3 At n r 



+ 



^® 'read 

g 2 n r At 2 



Note that the factor f drops out of the variance estimate. 

We now consider the case of n c co-additions. The variance then becomes, 



n z — ' 



Ae 

„2 



f.»4 ^ ^gTlj-At 2 



_ 1 ft (n 2 . - 1) 
3 At n r 



2<r 2 



g 2 n r At 2 



gn r n 2 At 2 
I 



! _ 1 ft (n 2 - 1) 



3 At n r 



gn c At 



1 1 St (n 2 1) 
3 At n r 



+ 



9_2 

, zcr read 

g 2 n r n c At 2 

2<X 2 . 

read 



g'n r n c 



At 2 ' 



(37) 
(38) 

(39) 
(40) 



(41) 
(42) 

(43) 
(44) 



Finally, let us define the effective read noise cr r ead,eff as 

V2&read 



&read, ef f 



\Jri r Ti c 



which has the units of electrons. The variance on I is then given by 



V, = 



gn c At 



1 ft (n 2 . - 1)' 
3 At n r 



+ 



read,ef f 

g 2 At 2 



(45) 



(46) 



Systematic effects produce additional terms in the equations for the variances. If these effects arc 
present, but not accounted for (as in our analysis), they will manifest themselves as an increased (i.e., larger 
than intrinsic) or perhaps variable "read noise" (see §4). 
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3. Continuous Sampling Technique 
3.1. Ideal Detector 

The notation given in §2 can also be used to estimate the pixel variances of an ideal detector read out 
with the "continuous sampling" (also known as the line-fitting or "up-the-ramp" ) technique (Chapman et 
al. 1990; Finger et al. 2000). In this technique, the array is read out repeatedly during the entire exposure. 
The count values from the readouts are then fit with a straight line and the slope gives the estimated source 
flux, I (in DN s _1 ). For n. r reads, each of duration 5t, during an exposure of At, the slope can be written as 
follows (Bevington 1969): 

nrE£it?-(E£iti) a ( } 

where tj is the time of the ith readout, and Sj is the count value recorded during this readout. This equation 
assumes that all readouts are given equal weight in the fitting process. If the readouts are equally spaced in 
time, tj = (i — f)St, At = (n r — l)5t, and the above equation can be simplified to yield, 

V™ r s (i - nr+1 ) 

I = 2 -'= 1 t{ 2_J (4g) 

a 

where a — n r (n r + 1) At/12 . (We refer the reader to the report by Sparks (1998) for the case where the 
readouts are not equally spaced in time.) For a perfectly linear detector exposed to a source with constant 
flux, the values recorded at the end of each individual read are given by 

Si = IU = (i- t)ISt . (49) 



3.2. Non-linearity corrections 

Non-linearity corrections are fairly straight-forward to implement with this readout technique. Again, 
we assume that the bias level is constant, the bias and dark current have been subtracted, and systematic 
effects have been accounted for before the non-linearity corrections are applied. Equation 49 yields the 
first estimate of the signal values, s\ for each read. The non-linearity corrected signal values can then be 
estimated in a manner similar to that for Fowler sampling, 

sf^sW-C^) • (50) 

Once all the signal readout values for all i reads have been corrected, the corrected slope can be re-determined 
from Equation 48. The procedure can then be repeated to yield more accurate values of s, and I, until a 
suitable convergence criterion is reached. 
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3.3. Variances 

From Equation (48) the variance in the slope recorded by a linear array can be calculated in a manner 
similar to that used to determine the uncertainties in MCDS: 



i=l 

= E 



91 V ~v^v- , J 9I \ ( 91 



(51) 



g g" 



»=i 

*EEf 

j=2 i=l S 



i- |K + 1) 



a 



+ 



i - i(n r + 1) 



i-i(n r + l) 



a 



(52) 



Using the definition of Sj above (Eq. 49), along with relations given in Equation 38 and the additional 
definitions 



we find 



E* 3 = 



n 2 (n + l) 2 



V 7 



' E z4 = ^( n+1 )( 2n + 1 )( 3n2+3n - 1 ) 



»=i 



read 



n r (n 2 — 1) 



120a 2 



6 I 



+ 



12a 2 

l2(T lead 



5 gn r At \n r + 1 / g 2 n r At 2 \n r + 1 

n r - 1 



6 I 



5 gn r At \n r + 1 



1 



6cr 2 



od,e// 



g 2 Ai 2 V n r + 1 



(53) 

(54) 
(55) 
(56) 



where we have used the definition of the effective read noise given earlier (Eq. 45). These equations are 
identical to that given by Garnett & Forrest (1993) once the different definitions of the exposure time are 
taken into account. However, they are substantially different from that derived by Chapman et al. (1990), 
who did not account for the correlation of the individual readouts in their analysis. 



4. Verification and Application 

To verify that the relation between observed count rate and variance predicted by Equation 44 for the 
MCDS technique yields realistic estimates of the variance in any given pixel, we obtained a series of flat 
field exposures with the SpeX instrument (whose array is read out using Fowler sampling) in the SXD mode 
(Rayner ct al. 2003). (Unfortunately, because we lack easy access to an array which is read out with the 
continuous readout technique, we were unable to obtain similar data to verify the procedures and equations 
given in §3.) Details about the readout procedures and electronics for SpeX can be found in Rayner et al. 
(2003) and references therein. Because SpeX is a cross-dispersed spectrograph, the flat fields exhibit a large 
range in observed count values, which makes them ideal for our purposes. Two series of exposures were taken 
with At = 17 s, St = 0.51 s, and successive n r values of 1,4,8,16,32,1 and 32,16,8,4,1,32. In all cases, 
n c = 1. For each combination of At and n r , nineteen flat field exposures were obtained, from which a mean 
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count rate value and a standard deviation were computed at each pixel. Non-linearity corrections were then 
applied to the individual frames and the means and standard deviations were re-computed at each pixel. 

For plotting purposes only, we computed the mean values and standard deviations of the variances in 
bins of the observed count values that were 50 DN wide. In Figure 3a, the mean variances are plotted versus 
the mean count values for one of the n r = 1 sets of flats. The results from the original data are shown 
as well as those obtained after applying linearity corrections. The same plot for one of the n r = 32 sets 
of flats is shown in Figure 3b. As can be seen from these figures, a linear trend of the variances of the 
non-linearity-corrected data with the observed counts, like that predicted by Equation 44, provides a very 
good representation of the results. 

Equation 44 was then fit to the actual (i.e., un-binned) flat field data, with a least squares procedure in 
which all data points were equally weighted, to determine the array-averaged gain g and read noise <7 reac i for 
each set of flats (i.e., for each set of At and n r values). We included only those pixels with values between 
< Stot/n r < 4000 DN in our analysis. Given the exposure parameters for the flat fields, this range of 
values represents most of the usable well for the SpcX array. Because the non-linearity corrections depend 
on n r , the corrections for the n r = 1 data sets are less than ~ 5%, while the non-linearity corrections for 
the n r = 32 data sets are less than about 15%, for S to t/ n r < 4000 DN. (These estimates have been derived 
from the SpeX non-linearity curve and simulations we have performed with SpeX flat field images.) SpeX 
documentation strongly advises users to keep values of S to t/n r below 4000 for all observations. The best fits 
are shown as solid lines on Figures 3a and 3b. Since most of the data points have values < 1500 DN, the 
fits are heavily biased to these lower values. (The distribution of Stot/n r values in the flat fields are shown 
at the top of the two plots.) 

The gain and readnoise values resulting from the fits are shown in Figure 4, plotted as a function of the 
number of reads. The values for the gain (Fig. 4a) were found to be reasonably consistent with one another, 
with a mean of 12.1 e~ DN -1 , although there is some indication of slight decrease in g as n r increases. 
We note that the smallest gain values were derived from the first images in both series of exposures, which 
were taken after the instrument and detector had been sitting idle for a long period. As seen in Figure 4b, 
however, there is a clear correlation between a rea d and n r . The variation in a rea d with n r is well described 
by a power law with an index of 0.16, 

a read w 36 • n° r 16 . (57) 

This relation implies that the effective read noise {<J r ead,eff, Equation 45) does not decrease as the square 
root of the number of reads, as might be expected, but rather decreases at a much slower rate, as n~ 0,34 . 
This is demonstrated in Figure 5 where we plot the effective read noise determined from the flat field data 
and the best fitting power law. 

An alternative description of the variation of the effective read noise with the number of reads is the 
following, 

&read,eff = A ■ n~ ' 5 + Constant , (58) 

where A is a constant proportional to the read noise <7 read - In this formulation, the read noise a rea d does 
have a constant value and the effective read noise does decrease as n~ , but there is an additional constant 
noise source that results in an additive offset. A fit of this form to the SpeX data is also shown in Figure 5, 
where it can be seen that it also provides an acceptable representation of the data. 

The variation of the effective read noise with n r determined from the SpeX data appears remarkably 
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similar to that presented in Figure 1 of Garnett & Forrest (1993). Although these authors claim that the 
effective read noise decreases as n~ 5 , close inspection of their figure reveals that a power law with n~ 0,5 
cannot reproduce the data points without an additive constant offset (Equation 58). A power law with a 
smaller index of about 0.4 would probably also fit their data well. 

As mentioned in §2, variation in the read noise with n r is a sign that systematic effects in the behavior 
of the array or the readout procedure are present but not accounted for in our analysis. Although there are 
several possible candidates for such an overlooked systematic effect (see §2), the exact cause for the observed 
variation of the read noise with the number of reads is uncertain at present. The variation is probably 
not a result of the non-linearity corrections, as the fits arc biased to a range of data values for which the 
non-linearity corrections are negligible. Furthermore, variations in the read noise with the number of reads 
have been observed before by others. For example, Finger et al. (2000) have suggested that the deviation 
from n~ - 5 in the effective read noise measurements of the Aladdin array installed in the ISAAC instrument 
at the VLT is due to additional Poisson noise from multiplexer glow and dark current. As the dark current is 
implicitly included in our formulation for the variance, increased dark current in the SpeX array is probably 
not the source of the additional noise. Inspection of dark frames generated by SpeX reveals no evidence of 
amplifier glow. Furthermore, it is straightforward to include an additional source term in Equations (5-9) to 
represent amplifier glow, which is present only during the readouts. This term leads to an additional noise 
term in the expression for the variance that increases as n r , which is significantly larger than the variation 
we observe. 

However, additional noise might result from the slight heating of the array due to the numerous reads. 
In order to reduce residual image effects the SpeX array is read out at 1 Hz when not integrating (Rayner 
et al. 2003). The act of reading out the array dissipates 20mW. An exposure consists of an initial burst (2 
Hz) of reads of the pedestal (heating), followed by a period of almost no power dissipation during the actual 
integration (cooling), terminated by a burst (2 Hz) of reads of the signal (heating). Although the array 
mount is controlled to 30.00 ± 0.01K, power is dissipated locally in the array, and so the variable heating 
effects of the readout might possibly cause much larger temperature differences. Changes in temperature 
(or the resulting instabilities) as a function of the number of reads would lead to a variation in the dark 
current rate and/or the bias level during the reads and therefore might explain the observed variation in read 
noise (see e.g., Boker et al. 2001). This may well be a common phenomenon in many readout techniques for 
near-infrared arrays. 

Because we have not explicitly included these systematic effects, or the corrections for non-linearity, in 
our equations, we should point out that the gain and read noise derived from our fits are not necessarily the 
"true", or intrinsic, values for the array. The non-linearity corrections of the individual images, for example, 
increase the variance of any given pixel value above that expected for the same pixel value recorded on an 
intrinsically linear array. Therefore, the estimated gain and read noise values determined from the fits will 
generally be smaller and larger, respectively, than the intrinsic values for the array. As far as data reduction 
is concerned, however, this does not pose a problem, as long as all images subject to the same systematic 
effects are reduced with the same parameter values. For example, the Spextool software first corrects each 
individual image for non-linearity before generating a variance image or proceeding with any reduction steps 
(see Cushing et al. 2003). The increased variances due to the non-linearity corrections are then correctly 
accounted for by adopting the gain and read noise values derived from our linear fits. Hence, these gain and 
read noise values are necessary for consistency in the data reduction routines. 



5. Summary 



We have derived fairly simple equations for the estimated variance as a function of the observed count 
rate at each pixel in a near-infrared array for correlated double sampling, multiple correlated double sam- 
pling, and equally-spaced continuous sampling read out techniques. We also give a general prescription for 
implementing corrections for non-linearity at each pixel. We have applied these equations to sets of flats 
obtained with the SpeX instrument at the NASA Infrared Telescope Facility and find that they provide a 
good representation to the non-linearity-corrected data. We also find that the read noise of the SpeX array 
varies with the number of reads as n° 16 , perhaps due to heating of the array as a result of multiple reads. 
Such heating may be generic to near-infrared arrays that employ (M)CDS as the readout technique. 
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Fig. 1. — Schematic of the signal recorded by an integrating node as a function of time on (a) a perfectly 
linear detector and (b) a real (i.e., non-linear) detector employing the MCDS read out technique. The two 
dotted vertical lines represent the starting time of the first read of the pedestal and the start of the first read 
of the signal. Each read is assumed to last 5t sec. Each solid vertical line represents the end of a read. The 
exposure time is At sec. 
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Fig. 2.— Spectra of Gl 846 from orders 6 and 7 of SpeX in the SXD mode (see Rayner et al. 2003). The 
spectra were reduced with the Spextool reduction package (Cushing et al. 2003). The spectra shown in the 
upper panel were reduced without applying the non-linearity corrections discussed in the text; the corrections 
were incorporated in the reduction of the spectra shown in the lower panel. Incorporation of the non-linearity 
corrections removes to a large degree the mismatch in both the flux levels and the slopes of the spectra in 
the overlap regions. 
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Fig. 3. — Measured mean variances and standard deviations as a function of count value for two sets of flat 
field data, (a) n r = 1 and (6) n r = 32. The values have been computed in bins of 50 DN wide. The open 
points are the original values, while the solid points represent the values measured after making non-linearity 
corrections to the individual frames (see text). The solid line represents a fit of Eq. 47 to the non-linearity- 
corrected data points. Since most of the data have count values below 1500 DN, the fit is heavily biased to 
match the variances in this range. The histogram of data values is shown at the top of both plots. 
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Fig. 4. — (top) Gain (in electrons/DN) and (bottom) read noise (in electrons) measured for the SpeX array 
as a function of the number of reads. The dashed line is the estimated mean value of the gain (12.1); the 
solid line is a fit of a power law to the read noise as a function of the number of reads. The power law index 
is 0.16. 
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Fig. 5. — The effective read noise as a function of the number of reads. The solid line is the best fit power 
law which decreases as n~ - 34 , significantly shallower than expected; the short dashed line is the expected 
decrease as n~ - 5 , normalized to the solid fitted curve at n r = 8. The long dashed line is the alternative 
fitting form in which the effective read noise decreases as n~ - 5 plus an additive constant. 



